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Abstract. We consider translation invariant quantum systems in thermo- 
dynamic limit. We argue that their energy-momentum spectra should have 
shapes consistent with effective models involving quasiparticles. Our main 
example is second quantized homogeneous interacting Fermi gas in a large cu- 
bic box with periodic boundary conditions, at zero temperature. We expect 
that its energy-momentum spectrum has a positive energy gap and a positive 
critical velocity. 
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1. Introduction 



1.1. Excitation spectrum of Fermi gas. In [8], one of the authors of this paper 
together with H. Cornean and P. Zih discussed a number of conjectures about the 
excitation spectrum of the interacting Bose gas at zero temperature with repulsive 
potentials in the thermodynamic limit. In particular, [8] conjectured that such 
systems have a quasiparticle-like excitation spectrum without an energy gap and 
with a positive critical velocity. These conjectures seem to be consistent with 
experimental data. In particular, they explain various phenomena related to the 
superfluidity. 

In this paper we would like to sketch a number of analogous conjectures about the 
interacting Fermi gas at zero temperature with an attractive interaction in therm- 
dynamic limit. We will argue that these systems should possess a quasiparticle-like 
excitation spectrum with a positive energy gap and a positive critical velocity. This 
conjecture implies in particular that the ground state energy is separated from the 
rest of the spectrum. In some situations, this ground state can be interpreted as 
a current carrying state, and plays an important role in the phenomenon of super- 
conductivity. 

The most robust quantity related to the excitation spectrum seems to be its 
infimum. Therefore, our main conjectures involve the infimum of the excitation 
spectrum (separately in the even and odd sector). They are based on the grand- 
canonical Hamiltonian for a fixed chemical potential. 

However, to have full information about the excitation spectrum, it is not enough 
to know its infimum. In fact, the HFB approximation suggests that the excitation 
spectrum of the Fermi gas has "lacunas" near its bottom. We make an attempt to 
express some conjectures about these lacunas. These conjectures are more compli- 
cated to state, probably also more delicate, and involve the canonical Hamiltonian, 
(that is, for a fixed number of particles). 

We do not prove our conjectures. However, we introduce a rather general class of 
model Hamiltonians for which these conjectures can be tested and, maybe, proven 
under some assumptions. These Hamiltonians consist of a two-body kinetic energy, 
not necessarily quadratic, and an interaction. The interaction does not have to be 
local (given by a local potential). We also suppose that our system has "internal 
degrees of freedom" (eg. "spin" ) . It seems important to assume that the interaction 
is in some sense "attractive", which means some kind of negative definiteness. This 
is suggested by the Hartree-Fock-Bogoliubov (HFB) method. 

We include in our paper short computations based on the HFB approach. The 
basic principle of these computations is well known, but in the literature they are 
usually presented in simple special cases. Our presentation applies to a rather gen- 
eral case. As a result of the HFB approach we obtain an approximate quasiparticle 
representation of our Hamiltonian with a dispersion relation possessing very special 
features: a positive energy gap and a positive critical velocity. We conjecture that 
this dispersion relation suggests basic qualitative features of the true excitation 
spectrum of interacting Fermi gas in thermodynamic limit. 
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1.2. Role of translation invariance. There exist many papers that study the 
energy spectrum of interacting Fermi and Bose systems. In particular, there ex- 
ist interesting works that study the HBF approximation in such systems. What 
makes our paper different is the role of translation invariance. This enables us to 
ask questions about the excitation spectrum, which we expect to have interesting 
properties. 

Many papers attempt to show that models based on quasiparticles give some 
kind of an approximation to realistic Hamiltonians, see eg. [9, 5]. However, only 
relatively crude features are considered in essentially all these papers. Typically, 
they study the energy or the free energy per volume in thermodynamic limit. We 
are interested in the excitation spectrum, which is a finer quantity and does not 
involve dividing by the large volume. The only rigorous result that we know devoted 
to the excitation spectrum of an interacting quantum gas is due to R. Seiringcr [21]. 
It concerns the Bose gas in finite volume and a mean field limit. 

We always assume that the interaction is translation invariant. In realistic phys- 
ical systems translation invariance is at most approximate. Superconducting mate- 
rials are perhaps the closest to idealized translation invariant models that we con- 
sider. Nevertheless, we believe that the picture presented in our paper is physically 
relevant also in many Fermi systems that are quite far from being translationally 
invariant, such as quantum dots and nuclei. In fact, one can argue that "traces" 
of translational degrees of freedom are present also in these systems, disguised as 
rotations and vibrations. In particular, the energy spectra of their odd and even 
sectors have various features consistent with our conjectures, see [3] for quantum 
dots and [19] Sect. 6.1 for nuclei. 

The immediate motivation of our paper is to state mathematically interesting 
rigorous conjectures together with heuristic arguments in their favor. Therefore, 
we do not strive at all costs to describe realistic concrete physical systems. The as- 
sumption of translation invariance helps to formulate a clean and rigorous definition 
of various concepts. 

Acknowledgement The research of J.D. and M.N was supported in part by the 
National Science Center (NCN) grant No. 2011/01/B/ST1/04929. The work of 
M.N. was also supported by the Foundation for Polish Science International PhD 
Projects Programme co-financed by the EU European Regional Development Fund. 

2. Quasiparticles and quasiparticle-like excitation spectrum 

The concept of "quasiparticle" , although often used, seems to have no satisfac- 
tory definition in the literature. In this section we attempt to give a number of 
rigorous interpretations of this term. We will also discuss spectral properties of 
quantum systems that can be described in terms of quasiparticles. 

The discussion of this section will be rather general and abstract. To a large 
extent it will be independent of the rest of the paper. 

2.1. Translation invariant quantum systems. The main object of interest 
of this paper are translation invariant quantum systems in thermodynamic limit. 
There are at least two approaches that can be used to describe such systems. 

In the first approach one starts with a construction of a system in finite volume, 
using A = [— L/2, L/2] d , the d-dimensional cubic box of side length L, as the 
configuration space. It is convenient, although somewhat unphysical, to impose the 
periodic boundary conditions, The system is described by its Hilbcrt space 7-L L , 
Hamiltonian H L and momentum P L . The spectrum of the momentum is discrete 
and coincides with ^Z d . After computing appropriate quantities (such as the 
infimum of the excitation spectrum) one tries to take the limit L — > oo. 
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Sometimes a different approach is possible. One can try to construct a Hilbert 
space W, a Hamiltonian H and a momentum P that describe the system on R d . 
This may be not easy, ft may require the use of refined techniques [7, 16]. It is 
probably not always possible. Note that in this case the spectrum of the momentum 
is expected to be absolutely continuous, with the exception of the ground state. 

The latter approach seems conceptually more elegant. Throughout most of this 
section we will adopt it. In most situations this will allow us to formulate some 
of the physical concepts in a concise manner. (Sometimes, however, it will lead to 
technical complications). 

In the next two sections we adopt the former approach, which is more down-to- 
earth. Thus only a family (H L , P L ) for finite L will be defined. 

To sum up, throughout most of this section by a translation invariant quantum 
system we will mean d + 1 commuting self-adjoint operators (H, Pi, ... , Pd) on a 
Hilbert space H. H has the interpretation of a Hamiltonian and P = (Pi, . . . , P d ) 
describes the momentum. 



2.2. Excitation spectrum. The joint spectrum of the operators (H,P) (which is 
a subset of will be denoted by sp(H, P) and called the energy-momentum 

spectrum of (H, P). 

We will often assume that H is bounded from below. If it is the case, we can 
define the ground state energy as E := inf spH. We will also often assume that H 
possesses translation invariant ground state <&, which is a unique joint eigenvector 
of H, P. In particular, H<P = E<5> and P$ = 0. 

Under these assumptions, by subtracting the ground state energy from the 
energy- momentum spectrum we obtain the excitation spectrum of (H,P), that is, 
sp(£f — E,P). We can also introduce the strict excitation spectrum as the joint 
spectrum of restriction of (H — E, P) to the orthogonal complement of $: 

Exc:=s P ((ff-i5,P)| w± ) (2.1) 

Thus if (E, 0) is an isolated simple eigenvalue of (H,P), then 

Exc = sp(H - E,P)\(0,0). 

Otherwise Exc = sp(H - E,P). 

We introduce also a special notation for the infimum of Exc: 

e(k) := inf{e : (e, k) G Exc}. 

The following two parameters have interesting physical implications. The first 
is the energy gap, defined as 

e := inf (sp(H - E) ) = inf{e(k) : k e R d }. 

Another quantity of physical interest is the critical velocity: 

c cr := mi -rz-r- 
Mo |k| 

Physical properties of a system are especially interesting if the energy gap e is 
strictly positive. In such a case, the ground state energy is separated from the rest 
of the energy spectrum, and hence the ground state is stable. 

Positive critical velocity is also very interesting. Physically, a positive critical 
velocity is closely related to the phenomenon of superfluidity, see eg. a discussion 
in [8]. 
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2.3. Essential excitation spectrum. One expects that most of a typical ex- 
citation spectrum is absolutely continuous wrt. the Lebesgue measure on R d+1 . 
However, it may also contain isolated shells continuously depending on the momen- 
tum. In this subsection we attempt to define the part of the excitation spectrum 
that corresponds to such a situation. 

Note that this is not easy in the abstract framework that we adopted in this 
section. Actually, in the next section, based on finite volume systems (H L , P L ), we 
will use a different approach to define isolated shells, see Subsect 3.7. 

We say that (e, k) £ belongs to Exc<j, called the discrete excitation spec- 

trum, if there exists S > such that the operator P has an absolutely continuous 
spectrum of uniformly finite multiplicity when restricted to 

Ran 1(\H - E - e\ < 5)1(\P - k| < S). 

The essential excitation spectrum is defined as Exc oss := Exc\ExCd- 

(We use an obvious notation for spectral projections of self-adjoint operators H 

and P: eg. t(\H — e\ < 5) denotes the spectral projection of H onto ]e — 8, e + S[). 
We introduce also a special notation for the bottom of Exc oss : 

e oss (k) := inf{e : (e, k) £ Exc oss }. 

Obviously, 

Exc D Exc css , 
e(k) < e oss (k), k£R d . 
Note that typically ExCd consists of a finite number of shells separated by lacunas. 
Abstract theory allows us to represent the Hilbert space T-L as the direct integral 
over R d given by the spectral decomposition of P, see eg. [6], 4.4.1. Suppose, 
in addition, that this direct integral can be taken with respect to the Lebesgue 
measure, so that we can write 

© 

H= J tf(k)dk. (2.2) 

keR d 

Then it is tempting to claim that 

sp(H-E,P) = |J sp(tf(k) - E) x {k}, (2.3) 

keK d 

Exc oss = ( |J sp css (H(k)-E) x{k}) Cl , (2.4) 

k£R d 

where sp ess denotes the essential spectrum and the subscript cl denots the closure. 
Unfortunately, at this level of generality there is a problem with (2.3) and (2.4). 
First of all, there is no guarantee that we can put the Lebesgue measure in (2.2). 
Secondly, the direct integral representation (2.2) is not defined uniquely, but only 
modulo sets of measure zero. 

In concrete situations, however, (such as quasiparticle systems considered in Sect. 
2.4) the direct integral (2.2) has an obvious distinguished realization involving the 
Lebesgue measure, for which the identities (2.3) and (2.4) are actually true. 

2.4. Quasiparticle quantum systems. Many important translation invariant 
quantum systems can be described in terms of quasiparticles, that is, indepen- 
dent bosonic or fermionic modes with appropriately chosen dispersion relations 
(the dependence of the quasiparticle energy on the momentum). 

Let us be more precise. For a Hilbert space Z, the notation T S (Z), resp. T a (Z) 
will stand for the bosonic, resp. fermionic Fock space with the one particle space 
Z. 
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By a quasiparticle quantum system we will mean (H{ T , Pf r ), where 

H h = Y, I ^(k)6*(k)6,(k)dk, (2.5) 

Pfr = Y [ H*(k)^(k)dk, (2.6) 

for some intervals Ii cR d , real continuous functions Jj 9 k >->• ojj(k), and creation, 
resp. annihilation operators 6*(k) and &i(k). Q is called the set of quasiparticle 
species and it is partitioned into Q s and <2 a - bosonic and fermionic quasiparticles. 

We are using the standard notation of the formalism of 2nd quantization: b* (k) 
and bi(k) satisfy the usual commutation/anticommutation relations. They are not 
true operators, only formal symbols, however the right hand sides of (2.5) and (2.6) 
are well defined as operators on the Fock space 

<8> r s (L 2 (/ 4 ))® ® r a (L 2 (/,)). (2.7) 

For i £ Q, the set Ii describes the allowed range of the momentum of a single ith 
quasiparticle and Wj(k) is its energy (dispersion relation) for momentum k £ R d . 
Note that Ii can be strictly smaller than M d - some quasiparticles may exist only for 
some momenta. This allows us more flexibility and is consistent with applications 
to condensed matter physics. It will be convenient to define 

/(k) :={i£ Q : kel,(k)} 

(the set of quasiparticles that may have momentum k £ R d ). 

Clearly, if we know the dispersion relations 3 k4 Wj(k), i £ Q, then we can 
determine the energy-momentum spectrum of (Hf T , Pf r ): 

sp(H fl ,P fl ) = {(0,0)} 

U {(wi 1 (ki) + "-+w in (k n ),ki + --- + k n ) : n = 1,2,3. ..} cl . 

Note that there is an obvious direct integral representation of the form (2.2) and 
the relations (2.3) and (2.4) hold. 

2.5. Properties of the excitation spectrum of quasiparticle systems. Let 

(H, P) be a quasiparticle system. The energy-momentum spectrum of such systems 
has special properties. First, we have 

(0,0) e sp{H,P), (2.8) 

because of the Fock vacuum state, which is a unique joint eigenstate of (H,P). 
Moreover, we have a remarkable addition property 

sp(H,P) = sp(H,P)+sp(H,P). (2.9) 

Assume now that the Hamiltonian (2.5) is bounded from below, or what is 
equivalent, assume that all the dispersion relations are non- negative. Then the 
Fock vacuum is a ground state satisfying E = 0, so that the excitation spectrum 
coincides with the energy-momentum spectrum. Thus we can rewrite (2.8) and 
(2.9) as 

(0,0) e sp(H-E,P), (2.10) 
sp(H-E,P) = sp(H -E,P)+sp(H -E,P). (2.11) 

Given (2.10), (2.11) is equivalent to 



Exc D Exc + Exc. 
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(2.12) 



Another remarkable property holds true if in addition the number of particle 
species is finite. We have then 

Exc oss = (Exc + Exc) cl . (2.13) 

Indeed, using the continuity of the momentum spectrum, we easily see that only 
1-particle states can belong to the disrete spectrum of the fiber Hamiltonians H(k). 

Before we proceed, let us introduce some terminology concerning real functions 
that will be useful in our study of quasiparticle-like spectra. Recall that a function 
R d 3 k e(k) is called subadditive if 

e(ki + k 2 )<e(k 1 ) + e(k 2 ). 

Let R d D I 3 k w (k) be a given function. Define 

f w (k) = inf{o;(ki) + ••• +w(k„) : ki + ••• + k„ = k, n = 1,2,3, .. .}, 
?es S ,«( k ) = inf{w(k!)H hw(k n ) : ki H h k„ = k, n = 2, 3, . . . }, 

(By definition, the infimum of an empty set is +oo). <r w is known under the name 
of the subadditive hull of u>. Equivalently, is the biggest subadditive function 
less than w. 

Note the relation 

?o,(k) = min{w(k) (k)}- 

Let us go back to a quasiparticle system (2.5), (2.6) with nonnegative dispersion 
relations. For k € R d , define 

^min(k) := mm{uj t : i G J(k)}. (2.14) 

Recall the functions e and e oss and the parameters e and c cr that we defined in 
Subsects 2.2 and 2.3. 

Theorem 2.1. (1) The bottom of the strict excitation spectrum is the subad- 
ditive hull ofui m i n : 

e(k) = ? Wmln (k), keM d . 

T/ie energy gap satisfies 

e = infw min (k). 

k 

(3) The critical velocity satisfies 

w min (k) . e css (k) 
c cr = mf — — — = mf — -— — . 
k#o |k| k#o |k| 

(4) If in addition the number of quasiparticle species is finite, then 

e css (k) = <r ess , Wmin (k), keM d . 

Note that we assume that the momentum space is R d . If we replace the momen- 
tum space R d with (that is, if we put our system on a torus of side length 
L) and we assume that all quasiparticles are bosonic, then all statements of this 
subsection generalize in an obvious way. However, because of the Pauli principle, 
not all of them generalize in the fcrmionic case. 
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2.6. Approximate versus exact quasiparticles. One often considers quantum 
systems of the form 

H = H b +V, P = Pfr, (2.15) 

where (H{ r ,P h ) is a quasiparticle system and the perturbation V is in some sense 
small. A description of physical systems in terms of approximate quasiparticles 
is very common in condensed matter physics. In particular, it appears naturally 
in the context of the so-called Hartree-Fock-Bogoliubov approximation, where one 
tries to optimize a quasiparticle description for a given quantum system [13]. 

Clearly, there is a considerable freedom in choosing the splitting of H into Hf T 
and V, and so quasiparticles of this kind are only vaguely determined. We will 
argue that in some cases a different concept of quasiparticles is useful, which is 
rigorous and in a way much more interesting. This concept is expressed in the 
following definition. 

Let (H, P) be a translation invariant system on a Hilbert space %. We will say 
that it is a quasiparticle-like system if it is unitarily equivalent to a quasiparticle 
system. 

2.7. Asymptotic quasiparticles. The above definition has one drawback. In 
practice we expect that the unitary equivalence mentioned in this definition is in 
some sense natural and constructed in the framework of scattering theory. 

Scattering theory is quite far from the main subject of this paper, which is mostly 
concerned with purely spectral questions. However, since it has been mentioned 
and is very closely related to the concept of a quasiparticle, let us give a brief 
discussion of this topic. 

For a number of many body systems the basic idea of scattering theory can be 
described as follows. Using the evolution e ltH for t —¥ ±oo, we define two isometric 
operators 

S ± : ® r s (L 2 (I t )) ® O r a (L 2 {Ij)) -> U. (2.16) 

are called the wave or M0ller operators and they satisfy 

HS ± = S ± H h , PS* = S ± P fl , 

where (Hf T , Pf r ) is a quasiparticle system. S := S + * S~ is then called the scattering 
operator. 

We will say that the system is asymptotically complete if the wave operators 
are unitary. Clearly, if a system is asymptotically complete, then it is quasiparticle- 
like. 

There are at least two classes of important physical system which possess a 
natural and rigorous scattering theory of this kind. 

The first class consists of the 2nd quantization of Schrodinger many body opera- 
tors with 2-body short range interactions [11]. One can show that these systems are 
asymptotically complete (see [10] and references therein). In this case the system 
is invariant wrt. the Galileian group and the dispersion relations have the form 
K a 3 k i-> E + ^— . Quasiparticles obtained in this context can be "elementary" 
- in applications to physics these are typically electrons and nuclei - as well as 
"composite" - atoms, ions, molecules, etc. 

Another important class of systems where the concept of asymptotic quasipar- 
ticles has a rigorous foundation belongs to (relativistic) quantum field teory, as 
axiomatized by the Haag-Kastler or Wightman axioms. If we assume the existence 
of discrete mass shells, the so-called Haag-Ruelle theory allows us to construct the 
wave operators, see eg. [17]. Note that in this case the system is covariant wrt. 
the Poincare group and the dispersion relation has the form R d 3kH \/m 2 + k 2 . 
Here, quasiparticles are the usual stable particles. 
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Let us stress that both classes of systems can be interacting in spite of the fact 
that they are equivalent to free quasiparticle systems. In particular, their scattering 
operator can be nontrivial. 

The above described classes of quantum systems are quite special. They are 
covariant wrt. rather large groups (Galilei or Poincare) and have quite special 
dispersion relations. 

2.8. Quasiparticles in condensed matter physics. The concept of a quasipar- 
ticle is useful also in other contexts, without the Galilei or Poincare covariance. 

An interesting system which admits a quasiparticle interpretation is the free 
Fermi gas with a positive chemical potential. We describe this system in Subsect. 
2.16. In this case the scattering theory is trivial: S + = S~ , and hence S = 1. 

It seems that condensed matter physicists apply successfully the concept of a 
quasiparticle also to various interacting translation invariant systems. 

One class of such systems seems to be the Bose gas with repulsive interactions 
at zero temperature and positive density. In this case, apparently, the system is 
typically well described by a free Bose gas of quasiparticles of (at least) two kinds: 
at low momenta we have phonons with an approximately linear dispersion relation, 
and at somewhat higher momenta we have rotons. This idea underlies the famous 
Bogoliubov approximation [4], see also [14, 8]. The phenomenon of superfluidity can 
be to a large extent explained within this picture. The model of free asymptotic 
phonons seems to work well in real experiments [18]. 

Another class of strongly interacting systems that seems to be successfully mod- 
elled by independent quasiparticles is the Fermi gas with attractive interactions 
at zero temperature and positive chemical potential. By using the Hartree-Fock- 
Bogoliubov (HFB) approach [19], which is closely related to the original Bardeen- 
Cooper-Schrieffer (BCS) approximation [1], one obtains a simple model that can 
be used to explain the superconductivity of the Fermi gas at very low temperatures. 
The corresponding quasiparticles are sometimes called partiholes. 

Note that the above two examples - the interacting Bose and Fermi gas - are 
neither Galilei nor Poincare covariant. This allows us to consider more general 
dispersion relations. However, we do not know whether these systems admit a 
quasiparticle interpretation or possess some kind of scattering theory. Unfortu- 
nately, rigorous results in this direction are rather modest. (There are attempts at 
scattering theory for some non-relativistic models of quantum field theory, see [12] 
and [15]. There exist also some results in a purely perturbative approach [20]). 

2.9. Quasiparticle-like excitation spectrum. The concept of a quasiparticlc- 
like system, as defined in Subsect. 2.6, is probably too strong for many applications. 
Let us propose a weaker property, which is more likely to be satisfied in various 
situations. 

Again, our starting point is a translation invariant system described by its Hamil- 
tonian and momentum (H, P). Let us assume that H is bounded from below, with 
E, as usual, denoting the ground state energy. We will say that the excitation 
spectrum of (H, P) is quasiparticle-like if it coincides with the excitation spectrum 
of a quasiparticle system (see (2.5) and (2.6)). Clearly, the excitation spectrum of a 
quasiparticle-like system with a bounded from below Hamiltonian is quasiparticle- 
like. However, a system may have a quasiparticle-like excitation spectrum without 
being a quasiparticle-like system. 

A quasiparticle-like excitation spectrum has special properties. In particular, it 
satisfies (2.10) and (2.11). 

There exists a heuristic, but, we believe, a relatively convincing general argument 
why realistic translation invariant quantum systems in thermodynamic limit at 
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zero temperature should satisfy (2.10) and (2.11). It was essentially described at 
length in [8], but for the convenience of the reader we reproduce it here. Note in 
particular that it the infinite size of the quantum system plays an important role 
in this argument. 

Consider a quantum gas in a box of a very large side length L, described by 
(H L ,P L ). For shortness, let us drop the superscript L. First of all, it seems 
reasonable to assume that the system possesses a translation invariant ground state, 
which we will denote by $, so that H$ = E$, P$ = 0. Thus (2.10) holds. 

Let (E+ei, kj) e sp(i/, P),i = 1,2. We can find eigenvectors with these eigenval- 
ues, that is, vectors <I>i satisfying H<&i = (E + e,)<I>i, P$i = ki$i. Let us make the 
assumption that it is possible to find operators Ai that are polynomials in creation 
and annihilation operator smeared with functions well localized in configuration 
space such that PAi ss Ai(P + kj), and which approximately create the vectors $j 
from the ground state, that is $i s=s A^. (Note that here a large size of L plays a 
role). By replacing $2 with e lyP $2 for some y and A2 with e lyP A2 e _lyP , we can 
make sure that the regions of localization of A\ and A2 are separated by a large 
distance. 

Now consider the vector <J>i 2 := AiAz^. Clearly, 

P$12 w (ki +k 2 )$12- 

$12 looks like the vector $i in the region of localization of Ai, elsewhere it looks 
like The Hamiltonian H involves only expressions of short range (the potential 
decays in space). Therefore, we expect that 

H$i2 w (E + ei + e 2 )$i2- 

If this is the case, it implies that (E + e x + e 2 ,k x + k 2 ) € sp(H,P). Thus (2.12) 
holds. 

2.10. Bottom of a quasiparticle-like excitation spectrum. Now suppose that 
(H , P) is an arbitrary translation invariant system with a bounded from below 
Hamiltonian. For simplicity, assume that its ground state energy is zero. Wc 
assume that we know its excitation spectrum sp(£f, P). There are two natural 
questions 

(1) Is sp(H, P) quasiparticle-like? 

(2) If it is the case, to what extent its dispersion relations are determined 
uniquely? 

In order to give partial answers to the above questions, recall the functions e and 
e oss , as well as the sets ExCd and Exc css that we defined in Subsects 2.2 and 2.3. 

Theorem 2.2. Suppose that the excitation spectrum of(H,P) is quasiparticle-like. 
Then the following is true: 

(1) e is subadditive. 

(2) We can partly reconstruct some of the dispersion relations: 

Exc d = {(Wi(k),k) : tee, keR d }\ Exc ess . (2.17) 
Consequently, for k satisfying e(k) < e css (k) ; 

e(k) = W m in(k), 

where w min was defined in (2.14). 

(3) If the number of quasiparticles species is finite, we can reconstruct e ess from 
e: 

e css (k) = inf{e(k 1 ) + e(k 2 ) : k = k!+k 2 }. (2.18) 

The existential part of the inverse problem has a partial solution: 
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Theorem 2.3. Suppose that M. d 3 k i-> w(k) fee a given subbadditive function. 
Consider the translation invariant system 



J w(k)fe k 6 k dk, P b = J k6 k 6 k dk. 



Hfr = 

Then 

e(k) = W (k), 
e oss (k) = inf{w(ki) + w(k 2 ) : k = ki+k 2 }. 

The answer to the uniqueness part of the inverse problem is negative. The only 
situation where we can identify dispersion relations from the spectral information 
involves ExCd, sec (2.17). The following example shows that we have quite a lot of 
freedom in choosing a dispersion relation giving a prescribed excitation spectrum. 
For instance, all the Hamiltonians below have the same excitation spectrum and 
essential excitation spectrum with e(k) = e css (k) = |k|: 

H= [ |k|(l + d|k| a )6 k 6 k dk, 

J\k\<c 

where c > 0, d > and a > arc arbitrary. 

2.11. Translation invariant systems with two superselection sectors. Sup- 
pose that a Hilbcrt space % has a decomposition W = V. + ® H~ , which can 
be treated as a superselection rule. This means that all observables decompose 
into direct sums. In particular, the Hamiltonian and momentum decompose as 
(H,P) = {H+,P+)® (H-,P-). Clearly, 

sp(H,P) = sp(H+,P + )Usp{H-,p-). (2.19) 

We will often assume that H is bounded from below and possesses a translation 
invariant ground state $ with energy E, which belongs to the sector % + . The 
sector H + will be called even. The other sector %~ will be called odd. 

Under these assumptions we will call sp(H + — E, P + ), resp. sp(£f~ — E, P~) the 
even, resp. odd excitation spectrum. We introduce also the strict even excitation 
spectrum: 

Exc+ := sp((ff + - E, P + ) ^ J (2.20) 

The strict odd excitation spectrum will coincide with the full odd excitation spec- 
trum: 

Exc" := sp(H~ -E,P~). (2.21) 

Finally, we define the even and odd essential excitation spectrum Exc^ s just as 
in Subsect. 2.3, except that we replace (H,P) with (H ± ,P ± ). 

We introduce also a special notation for the bottom of the sets Exc 1 * 1 and Excj^ s : 



Clearly, 

sp(H 



e ± (k) 


= inf{e : (e,k) G Exc 1 * 1 }, 




4 s (k) 


= inf{e : (e, k) G Exc^J. 




-E,P) = 


sp(H+ - E, P+) U sp(H~ - E, P~), 


(2.22) 


Exc = 


Exc + U Exc" , 


(2.23) 


Exc css — 


Exc+ ss U Exc ras , 


(2.24) 


e(k) = 


min{e-(k),e+(k)}, 


(2.25) 


e oss(k) = 


min Ks S ( k )>4 s ( k )}- 


(2.26) 
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2.12. Quasiparticle systems with the fermionic superselection rule. Con- 
sider a quasiparticle system (£ff r , Pf r ) on the Fock space (2.7). Define the fermionic 
number operator as 

iV a = ]T 6:(k)6i(k). 

Clearly, the fermionic parity (— 1)^ provides a natural superselection rule. If 
H = H + ® H~ denotes the corresponding direct sum decomposition, then the 
Hamiltonian and momentum decompose as 

(Hfr, P fr ) - (H+, P+) 8 (Hfr , Pfr ). (2.27) 

(2.27) will be called a two-sector quasiparticle system. 

Clearly, if we know the dispersion relations i, 9 k i->- w,(k), i g Q, then we can 
determine the even and odd energy momentum spectrum of (H^,P^): 

*p{H+,P+) = {(0,0)} 

U {(w il (ki) + ---+o;i B (k n ),ki + ■■■+!£„) : 

even number of fermions, n = 1, 2, 3 . . . } cl , 
s P (H fr ,P f r) = {(w il (ki) + ---+fa; jB (k n ),ki + "- + k n ) : 

odd number of fermions, n = 1,2,3... } c . 

2.13. Properties of the excitation spectrum of two-sector quasiparticle 
systems. Let (H, P) = (H + , P + )(B(H~ , P~) be a two-sector quasiparticle system. 
Clearly, we have 

(0,0) e sp(ff+,P+) (2.28) 

because of the Fock vacuum. Here are the properties of the even and odd excitation 
spectrum: 

sp(H + ,P + ) = sp{H+,P + )+sp(H + ,P + ) (2.29) 
D sp(H-,p-) + sp(H-,p-), (2.30) 
S p(if-,P-) = sp(if-,p-) + sp(tf+,P+). (2.31) 

Assume now that the Hamiltonian is bounded from below. Then the Fock vac- 
uum is a translation invariant ground state satisfying E = 0, so that the excitation 
spectrum coincides with the energy-momentum spectrum. Thus we can rewrite 

(2.28) -(2.31) as 

(0,0) e sp(H + -E,P+), (2.32) 
sp(H + -E,P+) = sp{H + -E,P+) +sp{H+ -E,P+) (2.33) 
D sp(H~ - E,p-)+sp(H~ - E,P~), (2.34) 
sp(H~ -E,P~) = sp(H~ - E,P-) +sp{H+ - E,P + ). (2.35) 

Given (2.32), (2.33)-(2.35) are equivalent to 

Exc+ D (Exc+ + Exc + ) U (Exc~ + Exc~), (2.36) 
Exc~ D Exc~+Exc+. (2.37) 

If in addition the number of particle species is finite, then 

Exc+ S = (Exc+ + Exc+) cl U (Exc - + Exc~) cl , (2.38) 

Exc ras = (Exc- +Exc+) cl . (2.39) 
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2.14. Two-sector quasiparticle-like spectrum. Consider now an arbitrary trans- 
lation invariant system with two superselection sectors {H,P) = (H + , P + )(B(H~, P~). 
We will assume that H is bounded from below and the ground state with energy 

E is translation invariant and belongs to the sector 

We will say that the excitation spectrum of (H + ,P + ) (H~ ,P~) is two-sector 
quasiparticle-like if it coincides with the excitation spectrum of a two-sector quasi- 
particle system. Such an excitation spectrum has special properties. In particular, 
it satisfies (2.32)-(2.35). 

There exists a heuristic general argument why realistic translation invariant 
quantum systems in thermodynamic limit should satisfy (2.32)-(2.35). It is an 
obvious modification of the argument given in Subscct. 2.9. 

Indeed, we need to notice what follows. (— 1) N * is always a superselection rule 
for realistic quantum system. In particular, if we assume that the ground state is 
nondegenerate, it has to be either bosonic or fermionic. We make an assumption 
that it is bosonic. 

The eigenvectors $i and $2, discussed in Subsect. 2.9, can be chosen to be 
purely bosonic or fermionic. Using the fact that the ground state is purely bosonic, 
we see that we can chose the operators A\ and Ai to be purely bosonic or fermionic. 
(That means, they either commute or anticommute with (— Consequently, 
we have the following possibilities: 

• Both $! and <E> 2 are bosonic. Then <J> 12 is bosonic. 

• Both $! and $2 ar e fermionic. Then <I> 1 2 is bosonic. 

• One of $1 and $2 is bosonic, the other is fermionic. Then $12 is fermionic. 
This implies (2.36) and (2.37). 

2.15. Bottom of a two-sector quasiparticle-like excitation spectrum. Sup- 
pose again that (H,P) = (H + ,P + ) © (H~,P~) is a translation invariant system 
with two superselection sectors. We assume that we know its excitation spec- 
trum. We would like to describe some criteria to verify whether it is two-sector 
quasiparticle-like. These criteria will involve the properties of the bottom of the 
even and odd excitation spectrum. 

Theorem 2.4. Suppose that the excitation spectrum of (H + ,P + ) © (H~,P~) is 
two-sector quasiparticle-like. 

(1) We have the following subadditivity properties: 

e"(ki+k 2 ) < e-(k 1 ) + e+(k 2 ), 
e + (ki+k 2 ) <e-(k 1 ) + e-(k 2 ), 
e + (ki+k 2 ) <e+(k 1 )+e+(k 2 ). 

(2) If the number of species of quasiparticles is finite, then we can reconstruct 
e- ss and e+ s from, e~ and e+ : 

e-Jk) = inf{e-(k 1 )+e+(k 2 ) : k = k!+k 2 }, 

e+ s (k) = inf{e+(k 1 ) + e+(k 2 ), e"(ki) + e"(k 2 ) : k = k!+k 2 }. 

2.16. Non-interacting Fermi gas. Let us give a brief discussion of the free Fermi 
gas with chemical potential \i in d dimensions. For simplicity, we will assume that 
particles have no internal degrees of freedom such as spin. 

The Hilbert space of n fermions equals T™ (L 2 (M d )) (antisymmetric square in- 
tegrable functions on (R d ) ra ). Let A^j denote the Laplacian A acting on the iih 
variable. Then the Hamiltonian equals 

n 

H n = £(-A (i) - /x). (2.40) 

i=l 
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It commutes with the momentum operator 

n 

pn : =E- iV w- 



i=i 



It is convenient to put together various n-particle sectors in a single Fock space 

T a (L 2 (A)) := K (L 2 (A)) . 
n— o 

Then the basic observables are the Hamiltonian, the total momentum and the 
number operator: 

H = 8 H n = [ a*(-A - ^)a x dx, 

n=0 J 

P = 8 P n = -i / a x V x a x dx, (2.41) 

n=0 J 
oo f 

N = © n = / a* a x dx, 

n=0 J 

where a* / a x are the usual fermionic creation/annihilation operators. 

The three operators in (2.41) describe only a finite number of particles in an 
infinite space. We would like to investigate homogeneous Fermi gas at a positive 
density in the thermodynamic limit. Following the accepted, although somewhat 
unphysical tradition, we first consider our system on A = [— L/2, L/2] d , the d- 
dimensional cubic box of side length L, with periodic boundary conditions. Note 
that the spectrum of the momentum becomes ^Z d . At the end we let L — > oo. 

It is convenient to pass to the momentum representation: 

H L = ^(k 2 - M )a*a k 

k 

P L = ^ka k a k , (2.42) 

k 
k 

where we used (2.41) and a x = L~ d/2 ^ k e ikx a k . We sum over k E ^Z d . 

It is natural to change the representation of canonical anticommutation relations 
and replace the usual fermionic creation/annihilation operators by new ones, which 
kill the ground state of the Hamiltonian: 

K '■ = a k> fe k := «k, k 2 > \i, 



Then, 



J k • — u k; 

6 k : = a k , 6 k := a k , k 2 < fi. 



H L = ^|k 2 -/i|fe k 6 k + J B L , 

k 

p l = ^>& k & k , 

k 

N L = ^sgn(k 2 -/x)6 k 6 k + C 7 



where 



e l = ]>> 2 - M ), 

k 2 <p. 

c L = E L 



k 2 <M 
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It is customary to drop the constants E L and C L . 

Set (temporarily) w(k) = |k 2 — fi\. In the case of an infinite space, the above 
analysis suggests that it is natural to postulate 



H = 



P 



N 



;(k)6^ k dk, 



k^k^kdk, 



sgn(k 2 - /i)&£6 k dk, 



(2.43) 
(2.44) 



(2.45) 



as the Hamiltonian, total momentum and number operator of the free Fermi gas 
from the beginning, instead of (2.41). 

The operators b^/b^ can be called quasiparticle creation/annihilation operators 
and the function k H> w(k) the quasiparticle dispersion relation. Thus a quasipar- 
ticle is a true particle above the Fermi level and a hole below the Fermi level. 

In Sect. 4 we describe a version of the BCS theory based on the Hartree-Fock- 
Bogoliubov approximation. This approximation suggests that the interacting Fermi 
gas can be described, at least approximately, by a Hamiltonian of the form (2.43) 
with a dispersion relation k h- > cj(k) that resembles |k 2 — except that its minimum 
is strictly positive. 

2.17. Examples of the energy-momentum spectrum. The energy-momentum 
spectrum of a Fermi gas described by (2.43) and (2.44) with various dispersion 
relations ui can sometimes have a curious shape. In the remaining part of this 
section we will illustrate this with several examples. We will present diagrams 
representing the energy-momentum spectrum. In the full and the odd cases, the 
dispersion relation w is a singular part of the spectrum and it will be denoted by 
a solid line. In the even case, the dispersion relation will be denoted by a dotted 
line. We will always consider the spherically symmetric case. 

First consider the non-interacting Fermi gas, which, as we argued above, has the 
dispersion relation w(k) = |k 2 — In dimension 1 its energy-momentum spectrum 
looks quite interesting: 




Figure 1. sp(H,P) in the non-interacting case, d = 1. 
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FIGURE 2. sp(H + , P + ) in the non-interacting case, d= 1. 




Figure 3. sp(H ,P ) in the non-interacting case, d = 1. 
Clearly, for d > 2 the energy-momentum spectrum is rather boring: 




FIGURE 4. sp(H,P), sp(H+,P+), sp{H~,P~) in the non- 
interacting case, d > 2. 
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In the case of an interacting Fermi gas, we assume that 

u = v/(k 2 - ) u) 2 +7 2 . (2.46) 

Calculations presented in Sect. 4, in particular equation (4.16), suggest that the 
dispersion relation obtained by the HFB method is qualitatively similar to (2.46). 



H 




*- P 



Figure 5. sp(_ff, P) in the interacting case, d = 1. 




1 ► p + 

Figure 6. sp(H + , P + ) in the interacting case, d = 1. 
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1 ► P~ 

FlGURE 7. sp(H~ , P~) in the interacting case, d = 1. 

H 




► |P| 

FIGURE 8. sp(H,P) in the interacting case, d > 2. 




i ► |P + I 

FIGURE 9. sp(_ff + ,P+) in the interacting case, d > 2. 
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Figure 10. sp(H ,P ) in the interacting case, d > 2. 



Again, the case d = 1 differs from d > 2. However, in all dimensions the energy 
gap and the critical velocity are strictly positive. 



In this section we describe some classes of interacting models that seem to be 
relevant for condensed matter physics. We also formulate conjectures about these 
models that we think are suggested by "physical intuition" . 

Internal degrees of freedom of particles, such as spin, play an important role in 
fermionic systems. They are in particular crucial in the BCS approach. Therefore, 
we will take them into account. We will assume that they are described by a 
finite dimensional Hilbert space C m . Thus the one-particle space of the system is 
L 2 (K d ,C m ). 

3.1. 1-particle energy. The kinetic energy of one particle including its chemical 
potential is given by a self-adjoint operator T on L 2 (R d , C™ 1 ). We use the following 
notation for its integral kernel: for $ s Is 3 (Mr, C m ), 



We assume that T is a self-adjoint and translation invariant one-body operator. 
Clearly, 

Ti lii2 (xi,x 2 ) = T i2)il (x 2 ,xi) 

= T il . i2 (x 1 +y,x 2 +y). 

The first identity expresses the hermiticity and the second the translation invariance 



We will sometimes assume that T is real, that is, invariant with respect to the 
complex conjugation. This means that Tj 1 ^ 2 (xi,X2) are real. An example of a real 
1-particle energy is 



3. The model and conjectures 




of T. 




where the ith "spin" has the mass rrii and the chemical potential /i^. 
If the operator T has the form 
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for some function t satisfying 



t(xi,x 2 ) = i(x 2 ,xi) 

= *(xi +y,x 2 +y), 

then we will say that T is spin-independent. 

Clearly, the 1-particle energy can be written as 

T M -( Xl ,x 2 ) = (27r)- d J TijW^-^dk. 

If it is real, then Tj ; j(k) = Tij(— k) 
If it is spin independent, then 

T i,i( k ) = T ( k )^.j- 
In the real spin-independent case we have r(k) = r(— k). 

3.2. Interaction. The interaction of the Fermi gas will be described by a 2-body 
operator V. It acts on the antisymmetric 2-particle space as 

(V$)i 1 , i2 (x 1 ,x 2 ) = ^2 / / V r 4l , i2 , l3iU (x 1 ,x 2 ,x 3 ,X4)$ J4i43 (x4,X3)dx3dx4, 

where $ £ Tl (L 2 (M d , C" 1 )) . We will assume that it is self-adjoint translation 
invariant. Its integral kernel satisfies 

Vii,i2,i3,i4 ( x l; x 2j x 3j x 4) = — ^» 2 ,»i ,13,14 ( x 2> x l i x 3 j x 4) 

= -^i,i 2 ,i4,i3( x l> x 2, x 4, x 3) 



(x 4 ,x 3 ,x 2 ,x 1 ) 

= ^i,i 2 ,i 3 ,u( x i +y,x 2 +y,x 3 +y,x 4 +y). 

The first two identities express the antisymmetry of the interaction, the third 
its hermiticity and the fourth - its translation invariance. We also assume that 
V(xi, x 2 , x 3 , X4) decays for large differences of its arguments sufficiently fast. 

We will sometimes assume that V are real, that means, they are invariant with 
respect to the complex conjugation. This means V^ li j 2i j 3j j 4 (xi, x 2 , x 3 , x 4 ) is real. 

We will say that the operator V is spin independent if there exists a function 
w(xi, x 2 , x 3 , x 4 ) such that 

ViiM, 13,14 ( x li x 2; x 3j x 4) 

= ^(«( x i, x 2,x 3 ,x 4 )(5 il ^ 4 (5 l2il3 - v(xi,x 2 ,x 4 ,x 3 )c5 ilA3 <5 J2 ,i 4 ), 



Note that 



u(xi,x 2 ,x 3 ,x 4 ) = w(x 2 ,xi,x 4 ,x 3 ) 



= w(x 4 ,x 3 ,x 2 ,xi) 

= u(xi +y,x 2 +y,x 3 +y,x 4 + y). 

It will be convenient to write the Fourier transform of V as follows 

y(xi,x 2 ,x 3 ,x 4 ) 

= (2tt)- M J c ikix 1+ ik 2 x 2 -ik 3 x 3 -ik 4 x 4 Q(k 4 , k 2 , k 3 , k 4 ) 

xS(ki + k 2 - k 3 - k 4 )dk 4 dk 2 dk 3 dk 4 , 

where Q(ki, k 2 , k 3 , k 4 ) is a function defined on the subspace k 4 + k 2 = k 3 + k 4 . 
(Thus we could drop, say, k 4 from its arguments; we do not do it for the sake of 
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the symmetry of formulas). Clearly, 

Q»i,»2,*3,*4(' £ i) ^2, k 3 , k 4 ) = — Qi 2) i 1 ,i 3 ,i 4 (k2, ki, k 3 , k 4 ) 

= -Qii,i 2 ,U,»3(kl,k2,k4,k 3 ) 

= Q*4,*3,»2,»i(k4)k3,k2,ki). 
If we assume that the interaction is real, then 

Qii,i2,i3,jfc(kl>k2,k 3 ,k 4 ) = Q»i,»2,i3,t*( — ki> — k 2 , — k 3 , — k 4 ). 
If we assume that the interaction is spin- independent, then 

<3iii 2 i3i4(ki,k2,k 3 ,k 4 ) = -(g(ki,k 2 ,k 3 ,k 4 )^ li4 ^ 2 j 3 - g(ki,k 2 ,k 4 ,k 3 )c) ili3 c5 i2i4 ), 

for some function q defined on ki + k 2 = k 3 + k 4 satisfying 

<7(ki,k 2 ,k 3 ,k 4 ) = g(k 2 ,k 1 ,k 4 ,k 3 ) 
= <?(k 4 ,k 3 ,k 2 ,ki). 

In the real spin-independent case we have in addition 

g(ki,k 2 ,k 3 ,k 4 ) = g(-ki, -k 2 , -k 3 , -k 4 ). 

For example, a 2-body potential V(x) such that V(x) = V(— x) corresponds to 
the real spin-independent interaction with 

u(x 1 ,x 2 ,x 3 ,x 4 ) = V(x! -x 2 )S(x 1 -x 4 )£(x 2 - x 3 ), 

9 (ki,k 2 ,k 3 ,k 4 ) = y'dqy(q)(5(k 1 -k 4 -q) ( 5(k 2 -k 3 + q). 

3.3. n-body Hamiltonian. The n-body Hamiltonian of the homogeneous Fermi 
gas acts on the Hilbert space T™ (L 2 (K d , C m )) (antisymmetric square integrable 
functions on with values in (C m )® n ). Let Tr^ denote the operator T acting 
on the ith variable and Vuj) denote the operator V acting on the (ij)th pair of 
variables. The full n-body Hamiltonian equals 

Hn = E T W+ E %)• 

l<i<n l<i<Cj <ti 

It commutes with the momentum operator 

n 

P n :=]T-iV Xl . 

i=l 

3.4. Putting system in a box. As discussed already in the previous section, to 
investigate homogeneous Fermi gas at positive density in thermodynamic limit it is 
convenient to put the system on a box A = [— L/2,L/2] d with periodic boundary 
conditions. This means in particular that the kinetic energy is replaced by 

T L (x u x 2 ) = ± J2 c ik '( Xl -^r(k), 
and the potential V is replaced by 

y L (xi,x 2 ,x 3 ,x 4 ) = 

k 1 ,... : k 4 e^Z d , ki+k 2 =k 3 +k 4 
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Note that V L is periodic with respect to the domain A, and V L (x) — » V(x) as 
L — > oo. The system on a torus is described by the Hamiltonian 

h l,u = J- T { % + J2 V (%) ( 3 - 2 ) 

l<z<n l<i<i<™ 

acting on the space T™ (L 2 (A, C m )) . 

3.5. Grand-canonical Hamiltonian of the Fermi gas. It is convenient to put 
all the n-particle spaces into a single Fock space 

r a (L 2 (A,C" 1 )) := © C(L 2 (A,C" 1 )) 

n=0 

with the Hamiltonian 



H L := © H 

n=0 



L,n 



= J a^T^ ^x^ - x, 2 )a X!l2 dxidx 2 

+ \ I j a xi,n a i2,i2 V i^i2^ 3 ,u( x i' x 2 J x 3 J x 4)ax3, 4 3 a x 4 ^ 4 dxidx 2 dx3dx4, 

where a Xji , a* ; are the usual fermionic annihilation and creation operators. The 
second quantized momentum and number operators are defined as 



P L := © P n ^ = -i / a* iV^ax.idx, 

n=0 



N L := © n = / at .-a„ ;dx. 

n=0 



Above we use the summation convention. In what follows we will usually omit the 
indices. 

In the momentum representation, 
H L = Yj T ( k K«k 

k 

+ 2^7 ^ Q(ki, k 2 ,k 3 ,k 4 )a ki a k2 a ka a k4 , (3.3) 

ki+k 2 =k 3 +k 4 

P L = ^ka k a k , 



k 



n l = $> k a k . 

k 

In the spin-independent case, the interaction equals 

9( k i' k 2 ) k 3,k 4 )a kiii a k2j a k 3 i: ,a k4;i 

ki+k 2 =k 3 +k 4 

In the case of a (local) potential, it is 

~^fd X! ^(q) a k+q,i a k'-qj a k'j a M- 
k,k',q 

H L ' ± (k) will denote the operator H L restricted to the subspace P L = k, (— 1) 
±1. 
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3.6. Infimum of the excitation spectrum. For a large class of potentials the 
finite volume Hamiltonians H L are bounded from below and have a discrete spec- 
trum. 

The ground state energy is defined as 

E L = inf spH L . (3.4) 

For k G we define the infimum of the excitation spectrum in the even/odd 

sector in finite volume: 

e L > + (k) := mfspH L '+(k)-E L , k ^ 0, 

6 L <+(0) := mf(sp(H L '+(0)-E L )\{0}), 

e L '-(k) := mispH L '-(k)-E L . 

For k € R d , we would like to define the infimum of the excitation spectrum in 
thermodynamic limit. To this end, first we define its finite volume version in a 
"window" given by 5 > 0: 

e L ^(k) := inf^^k'J : k' L e-Z d , |k-k' L |<<5}. 
Then we set 

e ± (k) = supfliminf (e L ' <5 ' ± (k))) . 

<5>0 ^ L^oo J 

Let us now formulate our conjectures about e . 

Conjecture 3.1. We expect that for a large class of potentials with attractive 
interactions the following statements hold true: 

(1) The functions M. d 3 k i-> e ± (k) <G K are continuous. 

(2) Let k e R d . Let (k s ,L s ) e f^Z d x [0,oo) obey k s k, L s oo. Then 
e ^±(k s )^e±(k). 

(3) If d > 2, then 

inf min(e~(k),e + (k)) = : e > 0. 

(4) If d > 2, then 

. min( £ -(k),6+(k)) 

mf ^ p- '- = : c cr > 0. 

k/o |k| 

(5) We have the following subadditivity properties: 

e"(ki+k 2 ) < e-(k 1 )+e+(k 2 ), 
e+(k 1+ k 2 ) <e-(k 1 )+e-(k 2 ), 
e + (ki+k 2 ) <e+(k 1 )+e+(k 2 ). 

To motivate the above conjecture, consider a model Hamiltonian 

H = J2 I Wi(k)&£ ii 6 M dk ) (3.5) 

where 6k, i, are (fermionic, but possibly also bosonic) annihillation/creation 
operators and /j 9 k4 Wj(k) are continuous functions defined on closed subsets 
Ii C M d . For k e R d , let cj m ; n (k) be the lowest dispersion relation defined as in 
(2.14). Assume that 

e := inf cj min (k) > 0, 

k 

. w min (k) 
k#o |k| 
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which is suggested by the HFB approximation, see Sect. 4. Then the infimuni of 
the even/odd excitation spectrum of the Hamiltonian H equals e ± (k) = <^ min (k) 
and has the properties described in Conjecture 3.1. 

Note that in this conjecture we expect statements (3) and (4) to be true only in 
d > 2. This is due to an argument based on the Galilean covariance in a box with 
periodic boundary conditions in one dimension explained in Sect. II B of [8]. It is 
valid both for bosons and fermions. 

3.7. Isolated quasiparticle shells. The quadratic part of the model Hamilton- 
ian (4.11) obtained as the result of the HFB approximation involves m fermionic 
quasiparticles (corresponding to the dimension of the "internal subspace" C m ). Its 
excitation spectrum will contain "lacunas" above its infimum separated by at most 
m shells. In this subsection we try to formulate an additional conjecture that takes 
these lacunas into account. This is more difficult than the conjectures from the 
previous subsection. It is also more dubious. 

Forke ^Z d andj>eN, we define the jth shell in finite volume in the n-body 
case 



i>j" n ' + (0) 



= the jth lowest eigenvalue of H L '"- + (k)-£ L , k ^ 0, 
= the j + 1st lowest eigenvalue of H L ' n ' + (0)-E L , 
= the jth lowest eigenvalue of H L ' n ~ (k)-E L . 



(Of course, when counting eigenvalues we take into account their multiplicity). 

Let N+ := {0, 2, 4, ... } and N~ := {1, 3, 5, . . . }. For k e R d , we would like to 
define the the jth shell in thermodynamic limit. To this end, first we define its 
finite volume version in a "window" given by 5 > 0: 

^' + (k) 

:= inf {i/""(k'A inf spH L n - E L < 5}. 

k^^Z",|k-ky<i,neR± 3 



Then we set 



Clearly, 



vf(k) = sup (liminf (vf' ' (k) 

5>0 ^ ^ 



^(k) = e ± (k), 
vf{k) < i^ +1 (k). 



Set 



e± s (k) :=sup{if (k) : j = l,2,...}. 
Let us now formulate the conjectures about ef ss . 

Conjecture 3.2. We expect that for a large class of attractive potentials the 
following statements hold true: 

(1) The functions R d 3 k i-> vf(k), e^ ss (k) € M+ arc continuous. 

(2) Let (k s ,L s ,n s ) e f^Z d x [0, ocfxN* obey k s k, inf spH L " n ° - E L ° ->• 0, 
L s — > co. Then 



^•"^(k 5 )^*f (k). 

A similar property holds for e^ s . 
(3) e^ s are related to ^ as follows: 

e" s (k) = inf{e-(k 1 ) + e+(k 2 ) : k = kx+k 2 }, 

e+ s (k) = inf{e+(k 1 )+e+(k 2 ), e" (k^ + e~ (k 2 ) : k = k 1+ k 2 }. 
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To justify this conjecture, let us note first that it is consistent with the spec- 
tral properties of the model Hamiltonian (3.5) if we assume that the number of 
quasiparticles is finite. 

We can try to be more precise. We expect that the functions Vj stabilize. In other 
words, for a certain m ± and j > m ± , all uf are equal to one another, and hence 
equal to ejf ss . Then it is natural to guess that the functions u)j that appear in the 
model Hamiltonian (3.5) and correspond to bosonic, resp. fermionic quasiparticles 
coincide with v- for j < m ± . 

Note that the HFB approximation, described in the next section, suggests that 
m + = and mT = to, where to is the number of internal degrees of freedom. In 
particular, this would mean that all quasiparticles are fermionic. This conjecture is 
probably too strong. One cannot exclude that the interaction leads to a formation 
of quasiparticles consisting of an even number of fermions. Such quasiparticles 
would be of course bosonic. 

4. The Hartree-Fock-Bogoliubov approximation applied to 
homogeneous fermi gas 

One can try to compute the excitation spectrum of the Fermi gas by approximate 
methods. Historically, the first computation of this sort is due to Bardeen-Cooper- 
Schrieffer. In its original version, the BCS method involved a replacement of qua- 
dratic fermionic operators with bosonic ones. We will use the approach based on a 
Bogoliubov rotation of fermionic variables, which is commonly called the Hartree- 
Fock-Bogoliubov method. Its main idea is to minimize the energy in the so-called 
fermionic Gaussian states - states obtained by a Bogoliubov rotation from the 
fermionic Fock vacuum. The minimizing state will define new creation/annihilation 
operators. We express the Hamiltonian in the new creation/annihilation operators 
and drop all higher order terms. This defines a new Hamiltonian, that we expect 
to give an approximate description of low energy part of the excitation spectrum. 

4.1. The rotated Hamiltonian. One can start the HFB method with a rotation 
of the fermionic creation/annihilation operators. For any k this corresponds to a 
substitution 



(* denotes the hermitian conjugation, # denotes the transposition and 7 denotes 
the complex conjugation). 

Here (4.2) guarantees that [a k ,ak] + = 1, (4.3) guarantees that [a k ,aL k ]+ = 0. 
Note that [a k , a k ']+ = 0, [a k , a*_ k ,] + = 0, k ^ k' are satisfied automatically. 

For a sequence 3 k # k with values in matrices on C m such that # k = #-k, 

set 



Note also that Ug is the general form of an even Bogoliubov transformation com- 
muting with P L . 



a k: = c k & k: + s kfr-k, a k = c k fek + s k frl k > 
where Ck and Sk are matrices on C m satisfying 

CkC k + s k s k = 1, 
c k sf k + s k c^ k = 0. 



(4.1) 



(4.2) 
(4.3) 




(4.4) 
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In this section we drop the superscript L, writing eg. H for H L . The Hamiltonian 
(3.3) after the substitution (4.1) and the Wick ordering equals 

H = B 

+ l -Y, °( k )^-k 0(k)6_ k 6 k + D(k)Kb* 

k k k 

+ terms higher order in b's. (4.5) 

Here are explicit formulas for B, D(k) and O(k): 

B = ^r(k)s k s k 

k 

+ 2l7 ^2 ^( k ' _k > ~ k ' k ') s kC- k C- k 'S k ' 

k,k' 

+ Jd ^2 < 5( k ' k '' k '' k ) s kSk'S k 'S k ; 



O(k) = 2r(k)c k s k 

+ Jd ^2 ^O*'' k ' _k ' k ) S k'C- k 'S_ k S k 
k' 




+ Jd X! *5( k ' k '' k '' k ) c k«k'Sk'Sk 



D(k) = r(k)c k c k - (r(k)s_ k .s_ k ) T 

+ J-[ ^2 < 5( k '' k '' k ' k ) S k'C- k 'S_ k C k 
k' 

^2 < 5( k ' k ' _k '' k ') c ks-kc-k's k / 

k' 

2 

+ £d X! *3( k ' k '' k '' k ) c kSk'S k 'C k 
k' 

2 ■> T 

--p 2^(<9(- k , k ',k', -k)s_ k s k /s k /s_ k ) . 

k' 

Note that the formulas for B, O(k) and D(k) are written in a special notation, 
whose aim is to avoid putting a big number of internal indices. The matrices c k and 
s k have two internal indices: right and left. We sum over the right internal indices, 
whenever we sum over the corresponding momenta. The left internal indices are 
contracted with the corresponding indices of r or Q. The superscript T stands for 
the transposition (swapping the indices). 



4.2. Minimization over Gaussian states. Let fi denote the vacuum vector, 
fie := Ugil is the general form of an even fermionic Gaussian vector of zero mo- 
mentum. Clearly, 

(n e \Hn e ) = b, (4.6) 

(btn \Hbtn e ) = B + D(k). (4.7) 
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Therefore, we obtain rigorous bounds 

E < B, 

E + er(k) < B + MD(k). 

We would like to find a fermionic Gaussian vector that minimizes B 
pectation value of H . We assume that there exists a stationary point ( 
B considered as a function of c and s. Bogoliubov transformations form 
hence the neighbourhood of the stationary point can be expressed in the 
way: 



Ck 


Sk 




Ck 


Sk 




"4 




Ck 






Ck 




.4 A . 



-the ex- 

•Sk,C k ) of 

a group, 
following 

(4.8) 



This means (including internal indices) that 



CiZ.k 

Ci;,k 
sum 



Cim,kC m /_k ^im,kS mZ,ki 
Sim,k'S m ; i k + ^im,]i c m l,kJ 
Cim,k s m l,k + Sjni.kC^ki 
^im,kC m / k Cim,k^ m/,k- 

We enter the above formulas into the expressions for B, 0(k) and D(k). 

We can always multiply Ck and Sk by a unitary matrix without changing the 
Gaussian state. Hence, we can assume that 

4 - ^-K)*A- (4-9) 

Since s' is a complex function we can treat s' and s' as independent variables. 
Ck = c k , Sk = Sk corresponds to s' = 0, s' = 0. Because of (4.9), we have 

d 



_d_ 



sl=0 
s'=0 



s>0 
s'=0 



0. 



0. 



Then, for example, taking the first term of B one gets 
d ^ 

/ /q|3,kW,k^a',k 



as: 



rt.k' 



S / =Q — 7"a / 3,k'C ar ,k'S | 8t,k', 
^=0 



which equals the first term of O(k) at c = c and s = s. Calculating other terms of 
B one finally gets 

sO(k)|c=g. (4.10) 



dB 
5s 7 



s'=0 



Thus the minimizing procedure is equivalent to O(k) = 0. This result is a special 
case of a more general fact discussed in [13] where it is called the Beliaev Theorem 
[2]- 

Thus, if we choose the Bogoliubov transformation according to the minimization 
procedure, the Hamiltonian equals 



H = B + D(k)6 k &k + terms higher order in fe's. 



(4.11) 



In the case of the model interaction considered by Bardccn-Coopcr-Schricffcr, 
described in many texts, eg. in [14], the minimization of B yields a dispersion 
relations that has a positive energy gap and a positive critical velocity uniformly 
as L — > oo, that is, 



inf DCk) > 0, 
k 



^ ^ > o- 

k#0 k 



(4.12) 
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This phenomenon is probably much more general. In particular, we expect that 
it is true for a large class of real, spin-independent and attractive interactions. In 
what follows we provide computations that seem to support this claim. 

Note that the reality and spin-independence of the interactions leads to a con- 
siderable computational simplification. By an attractive interaction we mean an 
interaction, which in some sense, described later on, is negative definite. 

Let us assume in addition that higher order terms in (4.11) are in some sense 
negligible. Then formally H is approximated by a quadratic Hamiltonian B + 
^ k Z?(k)6 k 6 k whose dispersion relation has a strictly positive energy gap and crit- 
ical velocity. We view this as an argument in favor of Conjectures 3.1 and 3.2. 



4.3. Reality condition. Let us first apply the assumption about the reality of the 
interaction. In this case, it is natural to assume that the trial vector is real as well. 
This means that we impose the conditions 

Ck = c_ k , «k = s-k- 
This allows us to simplify the formulas for B, D(k) and O(k): 

B = ^r(k)s k s k 

k 

+ ^fd X! ( 3( k ' k ' ~ k ' k ') s kC k C k <S k / 
k,k' 

+ Td^2 < 3^ k ' k '' k '' k )^kS k 'S k /S k , 
k,k' 

O(k) = 2r(k)c k s k 

+ Jd ^2 ^ k ' ~ k ' ~ k '' k ')( C kCk - S k S k )c k /S k / 
k' 

+ Jd ( 5( k ' k '' k '' k ) c k s k'Sk'S k , 
k' 

D(k) = r(k)(c k c k - s k s k ) 

+ Jj^2 _k ' ~ k '' k ') C kSkCk'«k' 

k' 

2 

+ ^Q(k,k',k',k)(c k s k /s k /c k - s k s k 'S k <s k ). 



4.4. Spin \ case. Assume that the "spin space" is C 2 and the Hamiltonian is spin 
independent. We make the BCS ansatz: 

a \ 1 " 

c k = cos^ k q 1 , 

.JO 1 " 
s k = sin0 k _ 1 , 
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where, keping in mind the reality condition, the parameters # k are real. Then 
B = ^V(k)(l-cos2(9 k ) 



^ ce(k, k') sin 26> k sin 26 k > 



4L d 



k,k' 



+ J]7 E ^ k < k ')(! - cos 20 k)(l - cos Wu), 



k,k' 



where 



Note that 



a(k,k') := -(<z(k,-k,-k',k')+<z(-k,k,-k',k')), 



/3(k,k') 



2<z(k,k',k',k)-g(k',k,k',k). 



Q (k,k') = a(k',k), /3(k,k') = /3(k',k). 
In particular, in the case of local potentials we have 



a(k,k') := -(y(k-k') + y(k + k')), 
/3(k,k') = 2t>(0) -y(k-k'). 



We further compute: 
O(k) 



D(k) 



(<S(k) cos2(9 fe + £(k) sin26» k ) 
(£(k) cos20 fe - <5(k) sin26» fc ) 



1 

-1 

1 
1 



where 



*( k ) = ^E a ( k > k ') sin2 ^ 

k' 

£(k) = r(k) + ^£/3(k,k')(l-cos20 k ,)- 



We are looking for a minimum of B. To this end, we first analyze critical points 
of B. We compute the derivative of B: 

d 29k B = 5(k) cos 2# k + f (k) sin 20 k . 

The condition die^B = 0, or equivalently O(k) = 0, has many solutions. We can 
have 



sin26» k = 0, cos20 k = ±l, 
or 



sin26» k = -e k , 7^ 0, 



cos 2# k = e k - 



e(k) 



(4.13) 

(4.14) 



V^ 2 (k)+^(k) 
where e k = ±1. 

In particular, there are many solutions with all # k satisfying (4.13). They cor- 
respond to Slater determinants and have a fixed number of particles. The solution 
of this kind that minimizes B is called the normal or Hartree-Fock solution. 

One expects that under some conditions the normal solution is not the global 
minimum of B. More precisely, one expects that a global minimum is reached by a 
configuration satisfying 

^ „„„o. ?( k ) 



sin 2# k 



v/^(k) + ^(k)' 



cos 20 k 



v /^(k)+e 2 (k)' 



(4.15) 
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where at least some of sin 26* k are different from 0. It is sometimes called a super- 
conducting solution. In such a case we get 

" 1 



D(k) = V£ 2 (k) + £2(k) 



1 



(4.16) 



Thus we obtain a positive dispersion relation. One can expect that it is strictly 
positive, since otherwise the two functions 8 and £ would have a coinciding zero, 
which seems unlikely. Thus we expect that the dispersion relation D(k) has a 
positive energy gap. 

If the interaction is small, then £(k) is close to r(k) and 8(k) is small. This 
implies that -D(k) is close to |r(k)|. If r(k) has a critical velocity for large k and 
D(k) has an energy gap, then this implies that D(k) also has a critical velocity. 

In other words, we expect that for a large class of interactions if the minimum 
of B is reached at a superconducting state, then D(k) satisfies (4.12). 

We will not study conditions guaranteeing that a superconducting solution min- 
imizes the energy in this paper. Let us only remark that such conditions involve 
some kind of negative dcfinitencss of the quadratic form a - this is what we vaguely 
indicated by saying that the interaction is attractive. Indeed, multiply the definition 
of 8(k) with sin20 k and sum it up over k. We then obtain 

sin 2 2fl kV ^ 2 (k) + e 2 (k) = --^ £ sin20 k a(k, k') sin 2^. (4.17) 
k k,k' 

The left hand side of (4.17) is positive. This means that the quadratic form given 
by the kernel a(k, k') has to be negative at least at the vector given by sin26'k- 

Let us also indicate why one expects that the solution corresponding to (4.15) 
is a minimum of B. We compute the second derivative: 

d 2 e k d 2 e k ,B = <J k , k , (- sin2M(k) + cos2# k £(k)) 
+ ^a(k, k') cos 26» k cos 20 k / 

+ ^/3(k,k')sin20 k sin20 k ,. (4.18) 
Substituting (4.15) to the first term on the right of (4.18) gives 

<5 k , k <V<5 2 (k) + £ 2 (k), 

which is positive definite. One can hope that the other two terms in the second 
derivative of B do not spoil its positive definiteness. 
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